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Improving Dynamical Domain- Wall Fermion Simulations * 

L. Levkova^ and R. Mawhinney^ 

'^Department of Physics, Columbia University, New York, NY, 10027 

We report on studies of the chiral properties of dynamical domain wall fermions combined with the DBW2 
gauge action for different gauge couplings and fermion masses. For quenched theories, the DBW2 action gives a 
residual chiral symmetry breaking much smaller than what was found with more traditional choices for the gauge 
action. Our goal is to investigate the possibilities which this and further improvements provide for the study of 
QCD thermodynamics and other simulations at stronger couplings. 
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1. INTRODUCTION 

For the lattice domain-wall fermion formula- 
tion there are two independent approaches which 
can be used to reduce the residual chiral symme- 
try breaking (mrcs) for practical sizes of Ls- One 
way is to modify the fermion part of the action 
and another way is to improve the gauge action 
used. In this paper we take the second approach 
by studying the DBW2 gauge action in dynamical 
DWF calculations and considering the effects of 
also adding an adjoint term to the gauge action. 

The DBW2 gauge action reduces TOres at finite 
Ls in quenched DWF simulations by suppressing 
the topological lattice dislocations in the gauge 
field configurations |ll2j . These simulations 
find that TOics is about two orders of magnitude 
smaller than the corresponding value for the Wil- 
son action at a given lattice scale. Here we inves- 
tigate the size of that improvement for the corre- 
sponding dynamical DWF cases at zero and finite 
temperatures. 

The form of the DBW2 action is as follows: 



tion where Ci — —1.4069 is nonperturbatively 
computed [21 • 

We extend the study of lattice dislocation sup- 
pression by adding an adjoint term to the DBW2 
gauge action for the dynamical DWF case in the 
hope that it will improve further mrcs in the 
strong coupling sector without affecting the long 
distance physics much. 

2. RESIDUAL MASS FOR DYNAMI- 
CAL ZERO-TEMPERATURE SIMU- 
LATIONS 

As a measure of the residual chiral symmetry 
breaking we use the ratio |2|: 

_ E..,(-^5,(2y,t)J5(^,o)) 



SdBW2 = 



iVc 



(1 -8ci) ReTrP^,(a:) 



Table 1 contains the parameters of the dynami- 
cal zero-temperature runs and the respective val- 
ues for TOros- The main feature of these results 
is that rrij-cs significantly increases with the gauge 
coupling. For the weak couphng of /3 = 0.80, 
mrcs is quite small and appears independent of 
the quark mass. For the strongest coupling from 
Table 1, /3 — 0.70, mres has risen to as much as 
one third of the input quark mass. 



+ Cl 



ReTr Rf,^{x) 



It is a renormalization group improved ac- 

*This work was conducted on the QCDSP machines at 
Columbia University and the RIKEN-BNL Research Cen- 
ter. LL and RM are supported by the US DOE. 



3. RESIDUAL MASS FOR DYNAMICAL 
FINITE TEMPERATURE SIMULA- 
TIONS 

For the finite temperature dynamical case the 
run parameters are shown in Table 2. We mea- 
sure mrcs from screening pseudoscalar correlators 



2 



rrif mp, x^/Dof m^, 



# conf. 



0.70 0.026 0.830(5), 1.4 0.0094(1) 32 

0.75 0.022 0.667(7), 0.8 0.00405(7) 41 

0.80 0.020 0.539(6), 1.6 0.00137(2) 75 

0.80 0.030 0.598(8), 0.8 0.00136(2) 72 

0.80 0.040 0.645(8), 1.7 0.00137(3) 32 

Table 1 

Parameters and mrcs of dynamical DWF zero- 
temperature calculations. All runs have volume 16"^ x 
32, Ls = 12 and Mg = 1.8. 

in the spatial z direction. Figure 1 shows the ratio 
function of the spatial coordinate z. Ide- 
ally after some z, mi-cs should be virtually a con- 
stant and the corresponding graph should have a 
plateau. In the case of A'^t = 4 and /3 = 0.75 the 
lack of a plateau in Wres (which rises noticeably 
with increasing z) likely indicates that the desired 
propagation in the four-dimensional boundaries 
is modified by propagation into the bulk fifth di- 
mension. This effect is noticeably reduced at the 
weaker coupHng of /3 = 0.80 and is the subject of 
further investigation. 



/3 TUf Nt nip ra-jr # conf. 

0.75 0.022 4 1.465(3) 1.344(4) 74 

0.80 0.020 4 1.472(3) 1.379(4) 37 

0.80 0.020 8 0.737(6) 0.53(1) 34 

Table 2 

Parameters of dynamical DWF finite temperature 
calculations. All runs have spatial volume 16^, Ls — 
12 and Ms = 1.8. 
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Figure 1. mrcs (2) for A''t = 4 and 8. 



crossings, which contribute to mrcs, are changing 
dramatically with /3. For the weak coupling of 
/3 = 0.80 the gap is open, there are few crossings 
and consequently mros is small. The situation for 
the stronger couplings worsens to the point that 
at the strongest coupHng of /? = 0.70 the gap 
is virtually closed and mrcs is about an order of 
magnitude larger. 

Another way to study the gauge configurations 
is to plot directly their plaquette distributions. 
Figure 3 shows the comparison of the plaque- 
tte distributions of the dynamical runs with the 
quenched DBW2 case. The conclusion is that the 
fermion determinant has the effect of pulling the 
distribution towards more negative plaquette val- 
ues, which are associated with lattice dislocations 
and larger mrcs- 
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From the results of section 2 and 3 it is clear 
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mrcs- We want to investigate the properties of the 

dynamical DBW2 gauge configurations which af- -0.2 ^,1,1,1 .""o 
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feet directly the residual chiral symmetry break- "z" - p — r ' 4^ ' r ^ 



ing. One way of doing that is by studying the 

spectral flow of the hermitian wilson dirac oper- "~ ^'^ * ®® wuuuuq^^ 

ator ^^Dw Figure 2 shows some typical spec- -0.2 P ,1,1,1 "Pq 
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tral flows for three of the zero-temperature dy- 
namical runs from Table 1. The two properties 

of the spectral flow, gap width and number of Figure 2. Spectral flows of 75 i^w for three dynamical 

zero-temperature runs from Table 1. 
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Figure 3. Plaquette distributions of dynamical runs 
from Table 1 compared with the quenched DBW2 
case. 



5. ADDING AN ADJOINT TERM TO 
DBW2 

Our idea to decrease mres in the strong coupling 
sector is to modify the DBW2 action by adding 
an adjoint term to it in the way shown below: 



5'dbw2 + adj — 



iVc 



(1 -8ci) 



ReTrPf^^a;) 



x;/j,<i/ 



+ ci ^ ReTri?^„(2;) 



where we choose Cf ~\- 2ca = 1, to keep the nor- 
malization for this three term action the same as 
for the plaquette action in the continuum limit. 
We expect that the modified action will suppress 
more strongly the negative tails of the plaquette 
distributions and thus it will get rid of the lattice 
dislocations which increase mios- 

Table 3 summarizes the results for TOrcs for 
three runs with the modified DBW2 action. The 
values for TOi-cs are an order of magnitude smaller 
than the values for the unmodified DBW2 with 
the same parameters. However to conclude that 
there is an improvement at a given lattice scale 
we need to measure the lattice spacings for the 
runs from Table 3. All of those runs turned out 
to be in the deconfined phase and to measure the 
scales correctly we need to repeat the calcula- 
tion at a larger volume. Nevertheless, Figure 4 
shows the effect of the adjoint term on the pla- 
quette distributions by comparing runs from Ta- 
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Figure 4. Plaquette distribution for two of the runs 
from Table 3 compared with the unmodified DBW2 
case with the closest in shape distribution. 

ble 3 with runs with the original DBW2 which 
have the closest in shape plaquette distribution, 
although the parameters are different. The modi- 
fied DBW2 renders the negative distribution tails 
shorter by comparison. If similar plaquette dis- 
tributions correspond to similar scales then there 
might be an improvement of roughly two times 
smaller m-ics for the action with the adjoint term. 
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ruf 


c/ 


Ca 




^ conf. 


0.65 


0.026 


1.2 


-0.1 


0.0031(1) 


168 


0.70 


0.022 


1.2 


-0.1 


0.00083(7) 


36 


0.75 


0.022 


1.2 


-0.1 


0.00045(3) 


59 



Table 3 

rriros results for dynamical DBW2 with adjoint term 
runs on 8^ x 16 volume, Ls — 12 and M5 — 1.8. 

6. CONCLUSIONS 

The DBW2 action significantly suppresses lat- 
tice dislocations and reduces mres at zero tem- 
perature for dynamical DWF runs at weak cou- 
plings. For N-t ~ A there is an unphysical increase 
in m,.cs- For Nt — 8, plots of mj-es are likely show- 
ing a plateau at weak coupling. The addition of 
an adjoint term to the DBW2 action in order to 
reduce further mres is still under investigation, 
but it shows promise for improvement. 
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